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Abstract 



We study the quantum dynamics of soliton-like domain walls in anisotropic 
spin-1/2 chains in the presence of magnetic fields. In the absence of fields, 
domain walls form a Bloch band of delocalized quantum states while a static 
field applied along the easy axis localizes them into Wannier wave packets 
and causes them to execute Bloch oscillations, i.e. the domain walls oscillate 
along the chain with a finite Bloch frequency and amplitude. In the presence 
of the field, the Bloch band, with a continuum of extended states, breaks up 
into the Wannier-Zeeman ladder — a discrete set of equally spaced energy lev- 
els. We calculate the dynamical structure factor S^^{q,uj) in the one-soliton 
sector at finite frequency, wave vector, and temperature, and find sharp peaks 
at frequencies which are integer multiples of the Bloch frequency. We further 
calculate the uniform magnetic susceptibility and find that it too exhibits 
peaks at the Bloch frequency. We identify several candidate materials where 
these Bloch oscillations should be observable, for example, via neutron scat- 
tering measurements. For the particular compound C0CI2 -21120 we estimate 
the Bloch amplitude to be on the order of a few lattice constants, and the 
Bloch frequency on the order of 100 GHz for magnetic fields in the Tesla range 
and at temperatures of about 18 Kelvin. 
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I. INTRODUCTION 



Bloch oscillationsETH of a quantum state in a periodic one-dimensional structure can be 
characterized as an oscillatory response to a constant force. The phenomenon is a remarkable 
example of the counter-intuitive nature of quantum mechanics. Not only does the particle 
oscillate in response to a static and homogeneous force, but the amplitude of the oscillation 
is inversely proportional to the magnitude of the force. 

The prototypical system in which to observe Bloch oscillations (bo) has historically 
been a band electron in an external electric field. In the absence of inelastic scattering 
and interband (Zener) transitions, the electron, in momentum space, continually traverses 
the Brillouin zone. Upon reaching a zone boundary, the electron is reflected to the oppo- 
site boundary and thereby reverses its momentum. In real space, the motion is likewise 
periodic — a constant force thus produces oscillatory motion. But the existence of BO has 
been controversial ever since its theoretical prediction many decades ago (Ref. ^ contains an 
early review). The controversy surrounding the theory was largely concerned with the use, 
or rather the misuse, of Bloch's theorem in systems where the lattice periodicity was explic- 
itly broken by a constant external field. Experiments were not conclusive because, in these 
early times, bulk solids were the only systems available. Their bandwidths, however, are 
typically too large, the lattice constants too small, and the inelastic scattering so frequent 
that the phase coherence of the particle state, essential for the existence of BO, is rapidly 
destroyed. 

However, with the advent of semiconductor heterostructures — in particular layered su- 
perlattices with large lattice constants, and thus small bandwidths, which allow many co- 
herent BO before phase coherence is lost — the debate about BO seems to be settled. Recent 
experiments now provide convincing evidence that BO of electrons can indeed occur. This 
confirmation comes from the detection of the coherent radiation emitted as ensembles of 
electrons execute BO.S There is also work which directly measures the physical displacement 
of electrons in superlattices as they oscillate.Ei 

The search for BO has not been limited to electrons. Recently, the effect was observed 
with ultra-cold atoms placed in an optical standing wave.i The standing wave served as the 
periodic potential, and a force was simulated by accelerating this periodic potential. The 
momentum distribution of the atoms had the time dependence expected from the theory of 

BO. 

Current studies of BO are moving beyond the establishment of their existence. In atom- 
optical systems, BO are being used to test other aspects of quantum theory. For example, 
nonexponential decay of unstable systems was very recently observed in the same systems 
used to observe BO and the Wannier-Stark ladder.0 In electronic systems, BO are being 
studied as a means of producing fast emitters of coherent radiation. The radiation emitted 
by the oscillating charges can be tuned over a wide range of frequencies simply by changing 
the electric field strength. Typical wavelengths are in the submillimeter range. 

An entirely different class of systems in which such coherence phenomena can be expected 
are magnetic systems, which brings us to the main subject of this work. It has been pointed 
out recently that BO should exist in purely magnetic systems — in the quantum dynamics 
of domain walls with soliton-like behavior.! This proposal for magnetic BO was based on a 
semiclassical treatment! of the quantum dynamics of extended domain walls moving in a 
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periodic potential and containing a large number of spins s, with s ^ 1. In the present work, 
we shall extend this investigation to the fully quantum mechanical regime of anisotropic 
spin-1/2 chains and demonstrate that BO occur in the quantum dynamics of elementary 
excitations such as spin solitons. Such solitons represent the extreme limit of a magnetic 
domain wall with a width of only one lattice constant. In contrast to the cases mentioned 
above (electrons and atoms), magnetic BO are an inherently many-body effect. The soliton 
motion is a cooperative phenomenon resulting from spin-spin interaction. Nevertheless, we 
will see that magnetic BO share many properties with their electronic counterpart. 

A remarkable feature of magnetic BO is that they give rise to oscillations of the magneti- 
zation at a Bloch frequency which can be continuously varied by an external magnetic field. 
Thus, besides being of fundamental interest, BO of magnetic solitons may also prove relevant 
for applications since they provide a natural source of magnetic dipole radiation — typically 
in the microwave regime. 

The outline of this work is as follows. In the next section, we define and discuss the 
one-soliton approximation in the presence of a magnetic field B applied along the easy axis, 
thereby extending the zero-field results obtained previously for antiferromagnets (AFMs)li2l 
and for ferromagnets (fms).0 We show that the external B field localizes the eigenstates 
of the Bloch band and discretizes the spectrum into a set of equally spaced levels which we 
call the Wannier-Zeeman ladder (wzl), in analogy to the Wannier-Stark ladder in electronic 
systems. We mention several spin-models capable of supporting BO, but our main focus is on 
Ising-like FM chains with biaxial anisotropy — a model which accurately describes compounds 
such aslli C0CI2 -21120. We then go on to discuss the spin correlation functions and calculate 
the wave vector and frequency dependent dynamical structure factor S^^{q, u) for two cases: 
zero and finite magnetic fields. For finite magnetic fields, the structure factor consists 
of peaks at frequencies corresponding to integer multiples of the Bloch frequency. These 
frequencies are typically in the GHz range. This result means that, for example, neutron 
scattering measurements on samples in thermal equilibrium can be used to observe the wzL. 
In zero magnetic field, we obtain a structure factor for FMs that is very similar in form to the 
one obtained for afms,0 apart from a wave-vector dependence that reflects the difference 
in magnetic ordering between FMs and afms. After a brief digression on the B —>■ limit, 
we turn to an investigation of specific materials which are promising candidates in which to 
observe BO of magnetic solitons. We focus on one particular Ising-like FM, CoCl2-2H20, 
which appears to be the best characterized of the ones we have identified. Other materials, 
though less well characterized, are potentially better candidates. Finally, we close with a 
summary of the main results, along with an outlook on future directions. 



In the following we concentrate on one-dimensional Ising-like magnets with nearest- 
neighbor exchange interactions J" and in the presence of a magnetic field b„ = g^s^n- The 
general anisotropic spin Hamiltonian is given by 



II. THE ONE-SOLITON APPROXIMATION 




(1) 
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where n denotes lattice sites and a denotes Cartesian coordinates. The exchange constants 
J" are either positive (fms) or negative (afms), with \J^\ ^ |J^'^|. Thus, the (Ising) 
easy-axis is along the z-axis. In systems with such a strong easy-axis, domain walls, or 
solitons (we shall use these terms interchangeably), are well defined. At sufficiently low 
temperatures, the system is in its ground state: ferromagnetic order for FMs, and Neel order 
for AFMs. The excitations consist of domain walls. In a pure Ising chain {J^ = = 0) 
with zero field {b'^ = 0), the spectrum consists of discrete energy levels, where each level 
corresponds to states with a fixed number of domain walls. If, as is usually the case, there 
are additional exchange couplings in the directions transverse to the Ising axis (J^ 7^ or 
7^ 0), then the degeneracy is lifted. The energy spectrum consists of a series of continuua 
separated by gaps. Each continuum consists of states with a fixed number of domain walls. 
The one-soliton approximation considers only the lowest band and neglects all transitions 
to higher bands. 

The spectrum described in the preceding paragraph has been verified numericallyEl for 
the anisotropic x-y model; for very large anisotropy (near the Ising limit), isolated bands 
with large gaps were observed, with the gap tending to zero as the isotropic limit was 
approached. The one-soliton approximation was used by Villain0 in his pioneering work on 
spin-1/2 solitons in Ising-like afm chains. The result of his calculation — the existence of a 
dispersive soliton mode (the Villain mode) below the two particle continuum — was further 
verified by theoretical^ and numerical^ work, inelastic neutron scattering experiments on 
the Ising-like AFMs CsCoCls (Refs. p!6| , p7D and CsCoBrs (Refs. p^ , pJD , as well as electron 
spin resonance,IH3 nuclear magnetic resonance ,0 and opticali experiments. 

Analytic work guided by physical reasoning, exact numerical work on finite systems, and 
experimental work on real physical systems all point to the existence of the Villain mode 
and thus justify the one-soliton approximation. What then are the conditions necessary 
for the Villain mode to exist? Two conditions should be met. First, a large easy-axis is 
required; the system must be near the Ising limit. This ensures not only that domain walls 
are well-defined excitations, but also that states with different numbers of domain walls are 
well separated in energy. Second, the domain walls must have dynamics. It is the dynamics 
of the domain walls which induce the band structure. In a pure Ising chain, for example, 
a localized domain wall is an eigenstate of the Hamiltonian and therefore has no dynamics. 
The Villain mode does not exist in pure Ising chains. 

Most of the existing work has focused so far on AFMs, and only recently has it been 
pointed out that dispersive soliton modes can also exist in Ising-like FMs0 and other spin 
chains, as discussed next. 



A. Primary Model 

In view of the candidate materials to be identified in Sec. 0, we focus mainly on spin-1/2 
FMs with a Hamiltonian given by (U) with J° > 0.0 We consider a static and homogeneous 
field b along the z-axis (the Ising direction) and rewrite the Hamiltonian in a more suggestive 
form: 

H = H^ + H'' + H^, (2a) 



4 



H' = --PT.S'nSn+i~b^S^^, (2b) 



= - jy) + s-s-^,\ (2c) 



= -\iJ' + J') E(^n + (2d) 

n 

where S*^ = 5^ ± iS^ are the usual raising and lowering operators. In preparation for the 
one-soliton approximation, we introduce the one-soliton states {\m,Q)}, defined by 

m m 

l^,i) = |---Tni---), |m,-i) = |---iiTT---)- (3) 

The right-hand sides are expressed in the S'^-basis. Here, m = 0, ±1,±2,... denotes the 
soliton position (m = corresponds to the center of the spin chain), and Q = ibl is the 
charge of the soliton. We can define the soliton position, charge, and translation operators — 
m, Q, and T„ respectively — as 

m\m,Q) = m\m,Q), Q\m,Q) = Q\m,Q), Tn\m, Q) = \m + n, Q) , (4) 

where m and Q are Hermitian operators, whereas T„ is unitary with = T^n- 

The one-soliton approximation is tantamount to considering our system as containing 
one domain wall and discarding those terms in (^ which create additional solitons. For 



example, H in (|2dD should be discarded because these terms will always create solitons, 
but H"' in ( p^ contains terms which translate the soliton by two sites, which will be identified 
with T2 and tJ^. Projecting (|]) onto the one-soliton subspace, we therefore obtain 

i^i-soi = ^J' + A(T2 + r_2)-&gm, A = {jy-r)/4, (5) 

where we have dropped the hats over the operators, and are measuring energy relative to 
the fully polarized (ferromagnetic) state. The band width A is different from zero only if 
Jx 7^ Jy, in contrast to afms or alternate field configurations (see below). In what follows, 
we shall work in a fixed charge sector. Thus, Q is effectively a constant and we will put 
Q = — 1 to be definite. Equation (||) is then formally equivalent to a single-band tight- 
binding model of an electron in an external electric field. Solitons play the role of electrons, 
and a magnetic field the role of the electric field. We can now go on to discuss the eigenstates 
and energy spectrum, both of which are qualitatively very different in the finite- and zero- 
field regimes respectively. Bloch oscillations can be derived either semiclassically from the 
zero-field solution, or fully quantum-mechanically from the finite-field solutions. We discuss 
each of these in turn. 
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1. Semiclassical Solution 



We will first consider the known physics in zero external field in order to make contact 
with previous work on AFMs,0 FMs,0 as well as with semiclassical derivations of BO. For 
6 = 0, and with periodic boundary conditions, the eigenstates of (|) describing the sohton 
are extended Bloch states labeled by the wave vector k: 

\k,Q) = -^Y.e''nm,Q). (6) 

V ^^tot m 

Here, we have set the lattice constant a eaual to one. The periodic dispersion relation 
resulting from these eigenstates is given bytH 

E{k) = ^r + 2Acos{2k). (7) 

This is the ferromagnetic analog of Villain's resultlll for AFMs. As mentioned below. Villain's 
result for the bandwidth A contains the sum, rather than the difference, of the transverse 
couplings. 

Semiclassically, we can reproduce the derivation of BO as it is given in conventional 
electronic treatments.iil In the absence of scattering, the effect of the field-dependent term 
bQm in is to drive the soliton through the band. The velocity v{k) of this motion is found 
by differentiating (0) with respect to k. On the other hand, the wave vector k acquires a 
time dependence through the force F = hk = b/a. Integrating v{k{t)) over time then yields 
the semiclassical Bloch oscillations. If x{t) denotes the soliton position, then 

x{t) = const. — cos(co'Bt), (8) 
with the Bloch amplitude Ab and the Bloch angular frequency ujb given by 

Ab = 4Aa/6, hujB = 2b = 2gfiBB. (9) 



2. Quantum Solution (Exact) 

We can compare the above semiclassical derivation of BO with a fully quantum treatment 
by keeping a nonzero magnetic field in the Hamiltonian right from the outset. For a spin-1/2 
chain with A^tot = 2A^+1 sites, and with 6 > 0, we can exactly diagonalize (|^) to yield energy 
eigenstates 

N 

\Em) = Cmn\n), (10a) 

n=-N 



Cmn = {n\Em) = J(„_„)/2(a), « = A/6, (10b) 



6 



where Jn is the ordinary Bessel function of order n, and we have dropped the label Q = — 1 
in the state vectors. The state is localized about the lattice site m in the sense that 
lim„^oo('^|-£'m) = for arbitrary m. We have thus chosen vanishing boundary conditions — a 
reasonable choice for localized states. The degree of localization is given by the argument of 
the Bessel function (see (|1^ below); strong fields act to localize the states. For example, if 
h —>■ oo, then a —>■ 0, and only the m = n term contributes to the sum in ([l0|) since J„(0) = 
5nfl- The eigenstate is therefore strongly localized. As the field decreases, the wavepacket 
spreads. The eigenstates at finite h are Wannier-like states and are thus qualitatively different 
from the extended Bloch states at zero field. This difference is also reflected in the energy 
spectrum: 

i^i-soi \Em) = E Cm,n [A {\n + 2) + \n - 2)) + bu \n)] 

n 

= E [A {Cm,n-2 + Cm,n+2) + hnCm,n] \n) 
n 

2 / 777, — Tl\ 

= ^~ ( — 2 — I ^^'^rn,n \n) = bmY^ Cm,n \n) = hm \Em) , (11) 



and thus 



= hm. (12) 



The important first term on the third line in (0) was obtained using the Bessel function 
identity@ zJy-i{z) + zJy^i{z) = 2i'J,y{z). The spectrum {Em} is discrete. It consists of a 
series of equally spaced levels with an energy-level spacing given by b. The analog of this 
spectrum for electronic systems is known as the Wannier-Stark ladder. Hence, for magnetic 
solitons, the term Wannier- Zeeman ladder (wzl) seems appropriate and we shall adopt this 
term below. 

The states {\Em)} are exact eigenstates only up to boundary terms: 

Hi.soi\Em) = bm\Em) + boundary terms, (13a) 
where a typical boundary term is given by 

(a) \±N) . (13b) 

2 

For large and finite b (hence finite a) , these boundary terms are negligible contributions. 
The soliton dynamics are not expected to depend on the boundary conditions for sufficiently 
large chains, especially since the soliton eigenstates \Em) are strongly localized by the field. 
For n = \m — N\/2 ^ 1, and fixed, finite a, the Bessel function in ( |13bD can be replaced by 
its asymptotic form@ 

J^(a)^-^(^y, (n»l), (14) 



V27m \2nJ 

and therefore decays as n"", provided the argument a remains fixed at some finite value. If 
N ^ 1, and assuming a is indeed fixed, then the correction terms in (|13D can be neglected 
so long as the center of the eigenstate wavepacket is not near the boundaries of the chain. 
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Exact numerical diagonalization of finite chains (see below) indicates that the states {\Em)} 
can be considered eigenstates everywhere except at a small, 6-dependent boundary layer on 
either end of the chain. The correction terms can hence be considered surface effects which 
become negligible in the thermodynamic limit. 

More severe problems arise in the limit of vanishing field (a ^ oo). In this limit, the 
above asymptotic form of the Bessel functions cannot be used since now both the order of the 
Bessel function (oc N) as well as its argument (oc 1/6^) diverge. The appropriate asymptotic 
form for J„(a) depends on the ratio n/a; for both a —>■ oo and n — > oo, we have@ 



Ua)^J-(-——A , , , J e-" (l«a<n). (15) 




This expression reduces to (|1^) in the limit a <^ n. Equation (|T^) shows that J„(a) is 
negligible only if \n'^ — a'^\ ^ 1. This effectively implies that no matter how large the 
system size becomes, there will always be a field sufficiently small, such that the present 
framework fails. In practice, this becomes problematic only when discussing the "Villain 
limit" of vanishing field, where, in any case, soliton collisions must explicitly be considered 
(see below). 

Another approach is to compare the analytic spectrum Em = bm with that obtained 
by numerically diagonalizing finite chains (using -?/i_soi)- In Fig- |l(a), we plot the results 
for b/A = 3. This shows that for such fields, our analytic expressions are quite good and 
can be used with confidence. By contrast. Fig. |I|(b) shows a comparison of numeric and 
analytic results for b/A = 5/N. We see here that the width of the boundary layer has 
increased greatly. (The boundary layer consists of those points which deviate substantially 
from the linear analytic result. In Fig. |l](a), the boundary layer is not discernible.) The 
boundary layer is not a function of A^. It is a function of b. As A^ increases with b fixed, 
the boundary layer therefore becomes less and less important. Nevertheless, as b tends 
to zero, the boundary layer increases until eventually the present framework of localized 
eigenfunctions must be abandoned. This is our first indication of the problems associated 
with the 6^0 limit. We shall return to this limit below in connection with the calculation 
of the dynamical structure factor. For most of this paper, however, we shall consider either 
sufficiently large fields such that Fig. |I](a) is the relevant scenario, or zero fields, where a 
dispersive mode with Bloch-like extended states is the correct description. These are the 
two experimentally relevant regions within the one-soliton approximation. 

Assuming a sufficiently large field (in the sense of the previous paragraph) and neglecting 
surface effects, the spectrum of (|]) is the WZL, Em = bm. The presence of the magnetic 
field thus destroys the continuous band structure of (0) and replaces it with an evenly- 
spaced ladder of energy levels, with the spacing between adjacent levels given by b. How 
are BO manifested within this fully quantum-mechanical framework? To compare with 
the semiclassical result of (^ and (H), we should specify an initial state and compute the 
expectation value of m as a function of time. Let us keep the initial state arbitrary and 
write \ip{0)) = Em^'ml^m)- Then, 

(m(t)) = {mnm) = e^ic™i' - IabRb (e— * Y.CmCn.^2) , (le) 
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FIG. 1. (a) Comparison of numerical and analytical results for the energy spectrum of @ for 
moderate field values. Only every fifth energy level is shown. (There is one level at each integer 
value of m.) For such fields, the analytic spectrum Em = bm is essentially exact, (b) The same 
plot as in (a), but with a much smaller field. Here, about 20% of the points show significant 
deviation between the analytic and numerical results. As h continues to decrease, this boundary 
region increases. For such low fields, the system exhibits approximate translational invariance and 
so a momentum representation becomes more appropriate. 



where Re denotes the real part. If we assume, for example, that all Cm are real, and also 
that X^m CmCm-2 = 1, then the oscillating piece of ([T6|) is identical to the oscillating piece 
of (P). However, if we take as the initial state one which is completely localized at one lattice 
site, |'?/'(0)) = |m), then the oscillating piece identically vanishes. This situation corresponds 
to the state vector evolving symmetrically about both sides of the initial position m in a 
sort of breather state. Thus, the behavior depends sensitively upon the initial conditions. 

Semiclassically, we have seen that BO result from the soliton being pushed unscattered 
through the band and undergoing Bragg reflection at the zone boundary. By contrast, in a 
full quantum treatment, BO result from the time evolution of a state which is not a (Wannier- 
Zeeman) eigenstate of the Hamiltonian. The frequency of oscillation is given by twice the 
energy between adjacent states. This is a result of the fact that (|^), with 6 = 0, contains 
an intrinsic periodicity of two lattice constants, which also accounts for the fact that the 
Brillouin zone is halved, as shown by (0). This effect can be understood semiclassically in 
terms of the Berry phase of the spins and is a result of the spin parity (see Ref. ^ for more 
details). 

Finally, we note that BO are a many-body effect which should be distinguished from 
Larmor precession of uncoupled spins in an external field. The former yields an oscillating 
magnetic moment along the direction of the applied field whereas the latter yields oscillations 
in directions transverse to the field. 

We have focused on biaxial FMs because the materials we have identified as candidates 
for observing BO are all biaxial FMs. In the next subsection, we show that BO of magnetic 
solitons can also exist in uniaxial FMs ( = J^) if the field is tilted away from the Ising-axis. 
We also show that BO exist in anisotropic AFMs if an inhomogeneous field is applied. 
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B. Other Models 



1. Ferromagnets 

It is not necessary to have a biaxial FM in order to create BO. In fact, a uniaxial FM 
[J^ = jy) may be preferable because the band width can then be externally controlled. For 
example, consider (0) with ferromagnetic couplings, with = J^, and with a homogeneous 
magnetic field along both the x- and ^-axes {z is still the Ising axis). The Hamiltonian is 
almost the same as (0); the differences are that if" now vanishes (because = J^) and 
there is a new term coming from the field in the x direction: 

H^' = -lb^j:iS:[ + S-). (17) 

This term gives hopping by one site. The one-soliton Hamiltonian is 

^r-soi = - h^T^ + T_i) - b^Qm, (18) 

which is practically identical to (^. Thus, uniaxial FMs exhibit the WZL if the external field 
is tilted away from the Ising axis. The important difference between this and the biaxial 
case is that now the strength of both terms are adjustable externally. The energy eigenvalues 



here are as in {^), but the eigenstates are replaced by 

\E^r' = Y.Jn.-n{-bVb^)\n). (19) 

n 

Although some of the materials we will identify in the following sections are reported to 
have only uniaxial anisotropy, the material which seems to the best characterized, and for 
which we provide the most detailed analysis, is one which is reported to be a FM with biaxial 
anisotropy. 



2. Antiferromagnets 

It is more difficult to achieve BO in afms because of the local Neel order. We mentioned 
above that BO can be viewed, at least semiclassically, as the result of applying a force on a 
particle in a band. How can one apply a force to an antiferromagnetic domain wall? The 
force is given by F = — VX^n^n ■ S„. If the external field b is homogeneous, then the force 
quickly averages to zero over the chain. However, applying an inhomogeneous field produces 
a net force. Equation (0) can still be used, but the couplings are now negative, and must 
be slightly altered to reflect the inhomogeneity of the field: 

n 

The one-soliton states must also be redefined. Rather than (^, we should write 

m 

\m,Q) = \---nninv--)- m 
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The charge Q can be defined by the first spin at the end of the chainJll Q = —1 for spin-up, 
and Q = +1 for spin-down. Now it is in (2c) which will always create additional solitons 
(and thus should be discarded in the one-soliton approximation). Conversely, H-^ translates 
solitons. As a specific example, we can take a magnetic field along the easy axis that linearly 
increases along the chain axis: 6^ = (such a field satisfies Maxwell's equation V ■ B = 0). 
Projecting down to the one-soliton sector, we obtain 

HilZ = ^1^1 + A^™(T2 + T_2) - ^6"Q(-1)" ("^ + ^) , = + J") /4. (22) 

This Hamiltonian is again similar to (^. The spectrum again consists of the WZL, but due 
to the antiferromagnetism, the dependence of the spectrum on the position variable m is 
slightly altered: 

Err. = \h\-ir{m + ]^. (23) 



The eigenstates are also very similar to (|10]). Only a must be slightly changed, again to 
refiect the antiferromagnetism: 

« = (-l)"^^^, (Q = -l). (24) 

Thus, much of what follows applies also to afms where the anisotropy can be either uniaxial 
or biaxial. One must only replace the homogeneous field with a linearly increasing field. 

In most of what follows, we shall consider biaxial FMs with static and homogeneous fields. 
The present subsection, however, shows that the same analysis can be carried over, almost 
without change, to uniaxial FMs and to biaxial and uniaxial afms. 



C. Conditions for Observation 

Here, we briefiy touch on the conditions necessary to observe BO in physical systems. 
This discussion will have to remain rather vague, as it requires a detailed knowledge of 
specific material properties and, based on that, further theoretical investigations. Still, we 
can list a few essential conditions in general terms, which are very similar to the ones studied 
in the context of mesoscopic effects in electronic systems.S 

First, there should be no Zener transitions (interband tunneling). The soliton oscillates 
only when it is refiected from one zone boundary to the opposite one within the same band. 
If the force on the soliton is too strong, it will gain so much energy at the top of the band 
that it will tunnel into a higher energy band. This tunneling will produce classical linear 
motion, rather than the quantum mechanical BO. Such transitions can be neglected if the 
Bloch frequency is much larger than the Zener transition rate. This effectively puts an upper 
bound on the field driving the particle. In the present context, it means that the exchange 
constant along the easy axis, a measure of the band gap, should be much larger than the 
magnetic field h. 

Second, there should be no inelastic scattering. Such events occur, for example, by 
emission and absorption of phonons, or via soliton-soliton interaction. Inelastic scattering 
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may destroy the phase coherent motion of the particles necessary for BO to occur. A detailed 
investigation into the nature and especially the magnitude of the spin-lattice coupling in the 
spin-chains we shall be discussing below is beyond the scope of the present work, and is, in 
any case, probably a matter best determined experimentally (e.g. from the measured line- 
width in the structure factor). It is known, however, that the spin- lattice couplings are far 
weaker than the analogous charge-lattice couplings. Also, if the soliton density is low enough, 
soliton-soliton interactions, being typically of short-range nature, can be neglected and the 
results of band theory are still valid. In Ising-like spin chains, a low-density requirement 
implies that the temperature should be less than the exchange coupling J^. Again this can 
be typically satisfied. 

Finally, elastic scattering, such as scattering from static random impurities, may also 
be a problem (although typically less restrictive). Here, one should consider the Anderson 
localization length induced by random disorder in low-dimensional systems. This length, 
which is on the order of the elastic mean free path of the propagating quasiparticle@ (in the 
present case, the soliton), should be greater than the Bloch amplitude, which places a lower 
bound on the B field driving the soliton. However, since the Bloch amplitude can typically 
be on the order of the lattice constant, this poses no severe constraints. 

Although the above conditions are demanding, it is quite encouraging that the pres- 
ence of extended states of dispersive sohtons fi.e. the Villain mode) has been established 
experimentally in Ising-like antiferromagnets.lii"ii This suggests that, at least for certain 
spin-chains, inelastic scattering and disorder can be neglected to first approximation. In 
the end, the inelastic mean free time Tin should be compared with the Bloch frequency ub- 
Bloch oscillations are possible if ujb > l/Tin. Typical values for ujb lie between 40-600 GHz 
(see below). 



III. THE DYNAMICAL STRUCTURE FACTOR 



In this section, we show that the dynamical structure factor S§^{q,uj) at finite field 
contains sharp peaks at integer multiples of the Bloch frequency ujb — clear evidence of the 
WZL. Thus, inelastic neutron scattering, for example, can detect the WZL. By contrast, 
we also calculate the dynamical structure factor SQ^{q,uj) at zero field and thus give the 
ferromagnetic analog of the Villain mode for afms. 

For a translationally invariant system such as the full Hamiltonian in (0), the dynamical 
structure factor is defined in the standard wayS 



S^\q,uj) = ^ r dte- 

ZTT J-oo 



■iujt 



{6SlM6Sm, 



(25) 



where 6S^ = — (S^). The Fourier transform of the spin operator is also defined in the 



■-q 

standard way: for a finite chain with A^tot = 2A^+ 1 sites, 5*^ = En=-N e*«"5^. (We continue 
to set a = 1.) If the eigenbasis is orthonormal and discrete, one may write 



MS'qli^m) 5{uJ-{En-E^)) 



s: 



(26) 



where is the energy eigenvalue. 

In the following subsection, we give results for 6 = and later give results for 6 7^ 0. 
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A. The Dispersive Mode (6 = 0) 



Without a magnetic field, the eigenstates are the Bloch states @. Substituting these 
states, for a fixed charge of Q = — 1, into (p6|) yields 



k,k' 



{k'\S'\k) 5{u - E{k') + E{k)) - 



s: 



(27) 



where E{k) is shown in (|^). (We neglect the constant factor of since it drops out in the 
end.) The partition function Z can be expressed in terms of the modified Bessel function of 
order zero: 



^ = Ee"^''^'^ = A^tot/o(2A/?). 



(28) 



The matrix elements of in the eigenbasis can be found by first noting that 



-(l/2)|m), m > n 
(l/2)|m), m < n, 



from which it follows that 



TV 



s:\k) = l\k)-^T.' 

V ^^tot m=n 



ikm 



m) . 



(29) 



(30) 



The one-soliton approximation breaks down as g ^ 0.11^1^3113 If the neutron transfers no 
momentum to the system, then the energy will likely go into the creation of another soliton 
(actually a soliton-antisoliton pair). These are not the processes we are interested in here 
and are not contained in our approximation. Rather, we wish to consider the process where 
a neutron scatters the soliton from one /c-state to another, which is less and less likely to 
happen as g — 0. Below, we shall give estimates of both the soliton density and coUision 
rate. 

Limiting the discussion to g 7^ 0, the matrix elements of are given by {k'\S^\k) = 
{G'"^Sk',k+q — e~*'^^5fcA;/)/(l — e*^), and so the modulus squared becomes 

2 1 



4sin^(g/2) 



{k'\ SI \k) 



From this, we can find the first term in (|27|) : 



{^k',k+q + ^kk') ■ 



(31) 



^o^(g,^) + 



s: 



(^(-) + I E - Eik + g) + Em) • (32) 



The Dirac delta function can be written as 

5{u - E{k + g) + E{k)) = 6{u + Qq sin(2/c + g)). 



(33) 



fin = 4Asing, 
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(34) 



where ^Iq represents an upper bound or cut-off on the frequency. Above this frequency, the 
structure factor vanishes. This is a direct result of the one-sohton approximation, which 
yields a dispersion relation with a finite bandwidth. For a given momentum transfer q, 
the maximum energy a soliton can absorb (or emit), while still remaining in the same 
band, is Qq. The delta function (^) can be written in a more usable form by using the 
relation 6{f{k)) = J2n^{^ ~ kn)/\f'{kn)\, where the fc„ are the zeroes of f{k). In our case, 
f'{k) = 2QqCos{2k + q) and the kn are fixed by the condition 

sin {2kn + q) = —uj/Qq. (35) 

If we substitute these results into (0), and in addition take the continuum limit {J2k 
Ntot Jdk /2tt and A^tot oo), we obtain 

S''(q,u)+ (S') 6(uj) = Tr^{S(uj) + ^y-^ — . (36) 

° \ ^/ ^ ' 4sm\q/2)\^ ' 27rZ V |2fig cos(2A;„ + g)| / ^ ' 

To perform the sum over n, we must first solve (P^f). Defining uo/VLq = sin0, where |0| < 7r/2 
since Ic^l < (|35|) is rewritten as sin(2/c„ + ?) = — sin0. The general solution is 

2kn + q={-lT^\n'K + <P), (37) 

which is valid for arbitrary integer n (n = 0, ±1, ±2, . . .). But not all values of n are allowed; 
the allowed values of n must be chosen such that < vr, implying that the allowed values 
of n depend on the values of both cj) and g, where \(j)\ < 7r/2 and \q\ < tt. For example, n = 
is always allowed, n = —3 is allowed only if vr < {4> — q) < 37r/2, and \n\ > 4 is never allowed. 
It turns out that for any (p and q, there are always four values of n which are allowed — two 
even values (either n = 0, 2 or = 0, —2) and two odd values (either n = 1,3, n = ±1, or 
n = —1, —3). Actually, the specific value of n makes no difference; the important point is 
whether n is even or odd (and there are always two of each allowed). For example, from ( P?]) 
we can write 



\2Qq cos(2fc„ + g) I = I {-lT2Qq cos (pl = 2jQl - u^. (38) 



Similarly, we have 



- 2A/3 cos(2A;„) = -2A/5 [(-1)" cose/) cos g - sin <^ sing] = | [±(-l)"+Yfi2 - cot g + 

(39) 

where the top sign (+) is for |g| < 7f/2, and the bottom sign (— ) for 7r/2 < |g| < vr. (These 
signs will also prove irrelevant.) We can use the above two results to write the sum in (|36D 

as 



Ttt^ r-, n= J =yexp -/3f±(-l)"+\/T]2^^cotg + cj 

^|2r],cos(2fc„ + g)| 2,/^]2_^2^ 2'" V ^ ^ V 



(40) 

We see now that the magnitude of n plays no role; it only matters whether n is even or 
odd. Since there are always two even and two odd values of n, the sum over the exponential 
becomes a hyperbolic cosine, and therefore the additional ± also becomes irrelevant: 



14 



Y^exp (±(-l)"+yr]2_^2cotg + cj)j = Ae'^'^^^osh (^pJ^^^lJ^ cotq) . (41) 



.2 

Gathering the above results, the structure factor becomes 

, ,2 I ( cosh fi/?./fi2 _cj2 cot q 

° ^ ^ \ ^ ^ 4sin2 g/2 I ^ ' 7^/o(2A/5)./^]f^ ' ^ ^ 



For the calculation of |(S'g)p, we proceed by noting that, for g 7^ 0, {k\Sq\k) = 
-6-^9^/(1-6^"). Since this is independent of fc, we have = | (fc|S'q^|A;) p = 1/4 sin2(g/2). 

Finally, subtracting this from (^2[), we obtain 5*0^(5, cj): 



5o^^(g,u;) = V^JL^ ^ (43) 

47rsin2(g/2)/o(2A/?)./fi2_cj2 



This expression is plotted in Fig. |^ for A = 0.925 K and T = 18 K. (These parameters 
are relevant for the following discussion on the material CoCl2-2H20.) The plot shows the 
frequency dependence at a fixed wave vector of 7r/2. In addition, we have convoluted the 
dynamical structure factor with a Gaussian R(uj) = [Ij ^/¥a')e~^ with ^Ja ~ 4 GHz. 
Thus, the square-root singularities at = ifi^ have been rounded, which is the expected 
effect of soliton-soliton collisions and other interactions we have not taken into account here. 
This structure factor is very similar to that found in AFMs by VillainS and Boucher et a/.,§ 
who also worked in the one-soliton subspace, and by Nagler tt a/. ,0^1 who worked in the 
two-soliton subspace. As argued by Nagler tt ai, similar results should be expected for all 
cases where the soliton number stays fixed. 

Equation (^31) assumes the existence of only one soliton in the system, whereas there will 
always be some finite density of solitons. Since the (thermal) energy required to create a 
soliton is J^/2, the result in ( ^3]) may be crudely weighted by a Boltzmann factor given byill 
SQ^{q,iu) — > 6-'^"^^/^S'o^(g, cij). A more proper treatment would be to include soliton interac- 
tion, with the possibility of creation and annihilation of solitons. Nevertheless, ( ^3]) should 
be qualitatively correct for Ising-like FMs just as the antiferromagnetic analog qualitatively 
describes the experimental findings. 

An important difference between the result above for FMs and the result for AFMs is the 
wave-vector dependence. The factor of sin^(g/2) in ( ^5D for FMs is replaced by cos^(g/2) 
for AFMs. This difference between sine and cosine is related to the difference in ordering 
between FMs and AFMs; for AFMs, g = vr is commensurate with the spatial spin order near 
the ground state, while for FMs, it is g = which is commensurate with the ordering. Thus, 
one can replace g by vr — g in going from FMs to AFMs; this changes sin^(g/2) into cos^(g/2). 
Actually, cot g and Qg also change signs, but (|i3| ) is invariant under this change. Also, the 
difference in the bandwidth A (J^ + for AFMs and — for FMs) has been discussed 
in Sec. fTWI. 



It should not be so surprising that the result for FMs and AFMs is so similar in the 
absence of magnetic fields. After all, we have already seen in Sec. that a variety of Ising- 
like models get mapped onto an effective tight-binding model for solitons in the one-soliton 
approximation. For example, in zero magnetic field, Eq. (^) for FMs is formally identical to 



15 



ST 




-25 25 
(B/2jt (GHz) 



FIG. 2. A plot of for q = tt/2, as given by (|4|), with A = 0.925 K and T = 18K. We 

have also convoluted SQ^{q,uj) with a Gaussian, and so the square-root singularities at w = =bi7g 
have been rounded, as expected if collisions and interactions are taken into account. 



Eq. (^) for AFMs. Differences between the two only arise in the presence of a field. But 
even then, one can choose different field configurations for FMs and AFMs in order to obtain 
similar structure factors; as shown before, the WZL also exists in AFMs if an inhomogeneous 
field is applied. 

Equation (|43| ) shows a divergence as g — *■ 0. This should not be viewed as a physical 
result, but rather as an indication of the failure of the one-soliton approximation in this 
limit. In contrast, when we consider the b ^ case below, we shall see that the g ^ limit 
is well behaved. This is because the soliton states are now localized and the one-soliton 
Hamiltonian is no longer translationally invariant on the lattice. (Thus, k is no longer a 
good quantum number.) This localization should dramatically decrease the collision rate. 
In this way, the magnetic field provides a physical cutoff for the above singularity at g = 0. 
We will come back to this issue in the following subsection. 



For an estimate of the soliton-soliton collision rate, we can employ the results in Ref. |T7. 
These authors have looked at the model described by -fff™! in (11)) with = 0, and so we can 
use their results for our ferromagnetic model if we simply substitute A for A^'"^. Following 
Ref. 0, the soliton density is given by rig = e"'^'^^/^ and the soliton occupation probability 
by p{k) = e~^^'^^^ /Z. The soliton velocity, given by the derivative of the dispersion relation, 
is Vk = — 4A sin(2/i;), and the average soliton velocity is defined as vq = (1/A^tot) J2k I'^'klpik), 
whose evaluation yields 



^^0 



4sinh(2A/3) 2A/3« 



^ SA/ir. 



7T/3Ioi2A/3) 

The collision rate uJc{k) depends on both and Vk, and is given by 



Uc{k) 



k' 



p{k') \v'j, -Vk\^ risVo 1 - (1 - tt/2) sm^{2k) 



(44) 



(45) 
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B. The Wannier-Zeeman Ladder (6 7^ 0) 



In this subsection, we derive a central result of this paper. In the presence of a magnetic 
field applied along the Ising axis, the dynamical structure factor will exhibit the signature 
of the WZL. We shall also find that the q ^ limit is well behaved, in contrast to the 
previous subsection on the zero-field regime. This, together with the fluctuation-dissipation 
theorem, will also enable us to calculate the uniform susceptibility x"{^)y which provides us 
with a measure of the magnetization autocorrelation function. We begin first with S^^{q,u!), 
followed by x"{^)- 



1. The Dynamical Structure Factor 
Using the eigenbasis of ( p!OD we can write the dynamical structure factor ( ]26|) as 



2N 



Sl'{q,u) = Z-^ ^ e-f'''^\{E^.^\Sl\E^.N)\ 5{uj - h{n - m)) - \{Sl)\ (46) 



m,n=0 



Here, we have shifted the origin of our coordinates to one end of the chain (with a similar 
shift in the partition function Z). Because of the simple energy-level structure ([T2|), the 
partition function is simply given by 



^ 1^^- ^^^^ 

We shall assume /?6A^tot ^ 1 in all that follows. 

To obtain the matrix elements of S^, note first that from equations ( pUj ) and (|T^), it 
follows that 

N 

(48) 



1 / ^ \ 

Sn\Em) = - ( \Em) — 2 ^ Cmm'\m') j . 

\ m'=n / 



Using this, along with the relations {Em\m') = Cmm' and {Em\Em') = Smm', we can find the 



matrix elements of S^: 



2 N N 



{Erh\Sg\Em) = -NtotSrhm^qO — ^ 6*^" X! ^mm'Cmm'- (49) 

n=—N m'=n 

The expansion coefficients Cmn are given in ( |10b| ). The second term on the right side can 
be brought into a more manageable form by interchanging the order of the summations: 

N N N m' 

n=—N m'=n m'=—N n=—N 

Performing the geometric sum over the exponential, the right side is rewritten as 
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1 ^ 

~ ^ m'=-N 



+JV)/2 

^ ^^^^^^ E (e-^''^-e^''^'"+^-^"^')) J-'(«)'^-'+('^-)/2(«)- (51) 

/ m' = {m-N)/2 



1 - 



For N ^ oo (but keeping m finite for now), tlie sum over tfie product of Bessel functions 
can be performed using the identity^ 



E Ur)Jk^n{p) \ '1^^^ {k^) = Jn{R) '^^^ \ M) (52) 



fc=— oo 



where, for all variables real and for n an integer, R and d are defined through the relations 
R = v^r^ + — 2rpcosip, R cos d = p — r cos ip, and R sin -d = r sin y). In particular, 

oo 

E e-2*^™'j„,(«)J^,+(^_„)/2(«) = (-2signg)('^-'")/2e«<?('^— )/2j^^_^^^^(2a|sing|). (53) 



m — — oo 



Now we can substitute (|5T|) and ( ^31) into (^), which enables us to write the matrix elements 
in closed form as 

Efh\S^\Eml = ( — ^tot^gO ~ ~^ ~ 1 ^fhm 

1 I / 1 \in—m Aqifh+m) 12 

^ (-?signg)^ -Y——J{.fh-m)/2{Q. (54) 

where C = (2A/6)| sing|. 

Taking the modulus squared and shifting the origin as in (^6|), we can evaluate the first 
term in the dynamical structure factor: 



/I 1 \ 

+ [^^-MO^T.^cos[q{m+l)]e-^'-j 



2sin^(g/2) 



The prime on the final summation indicates that only those terms with m — n even should 
be included. The first of three sums in (|^) is proportional to the derivative of the partition 
function: 

The second sum in (^) can be written as a geometric series by writing the cosine as an 
exponential. The result is 
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For the third and final sum in (|55|), progress can be made by breaking it up into one 
containing terms with n and m even, and another containing terms with n and m odd; for 
the odd sums, we put n — 2n + 1 and m 2m + 1, and for the even sums, n — 2n and 
m 2m. If we then recombine these two sums, and in addition take u = n — m, the final 
sum in ( |55D can be written as 



2N 



N 



N-m 



Y,' e-^'"^5{u - b{n - m)) ^ (l + e-^') J] e-^^'™ ^ 6{u - 2hv). 

m,n=0 m=0 v=—m 

Now we can interchange the order of the sums by using the identity 



(58) 



N N-m 
2f3bm 



m=0 



N N N-u 



5{uj - 2hv) = 5] 5(cu - 2hv) e 



-N 



Y: Siuj - 2bu) Y e- 

m=0 



u=l 



(59) 



The geometric sums can now be performed. Using also the fact that f3bN ^ 1, 
written as 



is 



2N 



e-'^^^^Ccu - 6(n - m)) = Ze^("-l"l)/2 Y S{uj-2bn). 

m,,n=0 



N 



n=-N 



Substituting 

5r(g,^) + 
1 



), and (|60|) into (^) yields 



+ 



^o(C) 



COS q — e 



-fib 



(l - e-^' 



(60) 



-iq p—f^b\ 



(3{uj-\uj\)/2 j2 /AX AT 



2sin^(g/2) 4sin^(g/2) 



N 



(61) 



The last term is the most interesting one; it will induce transitions between Wannier-Zeeman 
levels. Let us first, however, complete the calculation by determining |(S'g)p. From 
and (|5^), we have 



S: 



1 e~*^^ 

-iVtot^go - Yzr^ 



1 - e 



-13b 



-iq p—Pb 



+ 1 - 2Jo(C) 



cos q — e 



-(3b 



(1 



-(3b 



cosh(/36) - 1 



cosh(/36) — cosg / 4sin (g/2) 



• (62) 



Finally, from ( |6TD and (p^ ) we obtain the dynamical structure factor: 



N 



Sl\q.^) = 7^ Y Gn{q)5{uj-nujB) 

^ n=-N 



(63a) 
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FIG. 3. A plot of (H) for q = tt/2, with = a and Pujb = 0.4. For the material CoCl2-2H20, 
this corresponds to an applied field of = 0.81 T, u;B/27r = 154 GHz, and T = 18K. The peaks 
at a; = ztw^ are measured on the left vertical axis and the peak at w = is measured on the right. 
The peaks have been broadened by convoluting with a Gaussian as in (|6^), with y/a ^ 40 GHz. 



Gn 



cosh(/5ci;B/2) — cosg' 



(63b) 



Gn 



2sin^(g/2) 



n > 
n <0, 



(63c) 



where = 26 is the Bloch frequency. Equation (^) indicates in particular that inelastic 
neutron scattering is capable of mapping the Wannier-Zeeman ladder. A soliton initially in 
some given Wannier-Zeeman state may be excited to higher states. The neutron lineshape 
intensity for an excitation by n levels, is essentially given by square of the n-th order Bessel 
function J^(C)- This means that the relative amplitudes of the peaks can be controlled 
through the argument C, = (2A/6)| sing| by adjusting the external field b (hence also the 
Bloch frequency and amplitude). For example, in Figs. |^ and ^, we plot the structure factor 
as a function of for g = 7r/2. In these figures, we have convoluted the structure factor 
with a Gaussian: 



/oo 
dto' Riuj-uj')S'h\q,iu'), R{uj) 
-oo 



1 



(64) 



For the plot in Fig. ^ we have taken T = 18K, A = 0.925 K (both as in Fig. |^), and 
h = 3.71 K. The choice of these numbers is motivated by the candidate materials to be 
discussed in the following section. These values give a Bloch frequency of about 154 GHz and 
a Bloch amplitude of about one lattice constant.^ For the material C0CI2 -21120 discussed 
in the following section, this value of b corresponds to an external field of about 0.81 T, 
and the temperature of 18 K is just above the three dimensional ordering temperature. In 
Fig. 1^, we use the exact same parameters as in Fig. ^ except for b, for which we have put 
b = 3.71 K/4 = 0.9275 K. The Bloch frequency is correspondingly reduced by a factor of 
four, while the Bloch amplitude is increased by the same factor. The striking feature here is 
the relative amplitudes of the peaks, as compared with Fig. H with this smaller field, peaks 
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FIG. 4. The same plot as in Fig. |3|, but with b reduced by a factor of four. There is a similar 
decrease in the Bloch frequency, and a similar increase in the Bloch amplitude. The striking feature 
is the change in the relative amplitudes of the peaks, as compared with Fig. ^. Here we have plotted 
peaks up to three times the Bloch frequency all on the same scale. 



up to = ±3ujb can be distinguished on the same scale. The peaks in both figures away 
from uj = are the signature of the WZL. There exists one peak at every integer multiple of 
the Bloch frequency, with an amplitude given by the square of a Bessel function. 

We have repeated the above calculation numerically (for finite A^) and compared the 
results with the analytic ones just presented. For any given peak at u = uub {n = 
0, ±1, ±2, . . .), the numeric results converge to the above analytic ones as the system size 
2N + 1 grows. If instead we fix A^, then the numeric results converge to the analytic results 
as one moves away from the boundaries of the chain at ±iV. We may thus conclude that 
the numeric and analytic results agree in the thermodynamic limit. 



2. The Uniform Susceptibility 

It is interesting to note that in contrast to the dispersive mode for 6 = 0, the q ^ limit 
is well behaved when 6 7^ 0: 

We have verified this result by performing the calculation at g = from the outset, as 
well as by calculating the imaginary part of the zero wave- vector susceptibility x"{^) the 
Matsubara formalism, and then using the fluctuation-dissipation theorem: 

X"{^) = ^^(1 - e-^nS','{u;) = ^ (1 - e-'--) [S{uj - ub) - S{u + ub)] ■ (66) 

Here, we have expressed x"{^) iii units of /x^ x sec. Equation ( |66D represents the response 
from only one soliton. There should be one such factor for each soliton in the system. 
Assume, for example, that one soliton exists every ten lattice sites for each chain in the 
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sample. (This requires the Bloch amphtude to be less than ten lattice sites.) A single- 
crystal of CoCl2-2H20, with a volume of Imm^ will then contain up to 10^® solitons. Thus, 
the signal can be quite large and should thus be observable in standard magnetization 
measurements (using, for instance, cantilever or SQUID technology). 

Again, the structure factor and the susceptibility should be observable as long as the 
inelastic scattering rate is less than ub- Concrete estimates for the soliton collision rate, 
for example, in the presence of a 6-field are even more difficult to obtain than for 6 = 0. 
Nevertheless, the rate for sufficiently large b should be far lower than the rate for 6 = 0. 
Indeed, when we turn on the field, the solitons become localized and execute BO about their 
mean positions. By tuning the field appropriately, the Bloch amplitude can be kept much 
smaller than the average distance (i.e. inverse density) between the localized solitons, and 
in this case the soliton-soliton interaction can be expected to be negligible. 

We can now go on to compare the zero-field limit of the above "Wannier-Zeeman" results 
with the previous "dispersive" results. 



As the field decreases, the Bloch amplitude increases. At some point, the Bloch amplitude 
will be equal to the spacing between solitons in the chain. At this point, collision effects, 
which are not included in our theory, become important and BO will be suppressed. This 
means that, for a meaningful comparison between our approximate theory and experiment, 
we must consider either the zero field regime {B = 0), or the regime where B is sufficiently 
large such that our one-soliton approximation is justified. 

In this sense, the limit of i? — > within the one-soliton approximation, has no experi- 
mental significance. However, it is still interesting from a technical point of view since this 
limit shows some features reminiscent of the classical limit of a quantum system, in the 
sense that there might be no pointwise convergence.^ Other well-known examples are the 
harmonic oscillator and a particle in a linear potential.0 In general — and especially when 
interference effects are important, as in the BO problem — some, usually ad hoc, averaging 
procedure must be employed in order to obtain a meaningful classical limit. 

The problems show up already at the level of the eigenfunctions and eigenvalues. At 
finite field, the eigenfunctions are the localized Wannier states ( |10bD , {n\Em) = [1 + 
(— l)"*"*^)] J(m-n)/2(tt)/2, which satisfy vanishing boundary conditions. But at 6 = 0, the 
extended Bloch states {n\k) = e'^^'^/\/N^ satisfy periodic boundary conditions. It is thus 
not too surprising that as 6 — > 0, the Wannier states do not converge pointwise to the Bloch 
states. The same holds for the energy eigenvalues; the band structure cannot be recovered 
from the WZL in the 6 — limit. 

Similar remarks apply to the partition functions. Keeping a finite system size A^tot = 
2A^ + 1, the partition functions for finite and zero fields are respectively given by 



C. Discussion of Results (b 0) 



N 





sinh(/36/2) 



(67) 



m=-N 



IT 



2A/3 cos(2fc) 



Artot/o(2A/5) . 



(68) 



k=—7r 
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These expressions imply that the 6 — hmit must be coupled to the A^tot oo limit. For 
example, if we are in the regime where A/? -C 1 and pb -C PbNtot "C 1, then to second order 
in both Ap and pb, Z^, converges to provided we make the identification bNtot = 2a/6A. 
On the other hand, if we are in the regime A/5 <^ 1 and /96 <^ 1 <C PbNtot, the partition 
functions cannot be matched. 

Next, we look at the dynamical structure factor (|63|) which we rewrite here in a slightly 
different form: 
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0, the second term tends to zero and so we need only concern ourselves with the first 



The Bessel function J|^|/2b (C) 
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J|^l/2b (l^gl/26) can be expanded in its asymptotic 
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It is the rapidly oscillating factor which prevents the pointwise convergence in the 6 — > 
limit; some averaging procedure is required. To illustrate this, let us crudely replace the 
oscillatory factor with 1/2 (which is accurate only for u;^ <^ Q"^): 
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In the 6 — > limit, the sum over n in (|69|) can be made continuous 
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Inserting ( fflf ) and (|72|) into (^), and assuming |ct;| < \^q\, we can write 



''-o 47rsin2(g/2)./fi2-cu2 



(73) 



Although the last two expressions are not equal, they are nevertheless quite similar. In 
Fig. H, we plot the two results, Eqs. (^) and (^3]), using the same parameters as in the 
previous figures. The small deviations can be traced back to the rapidly oscillating sine 
squared function in ([7D|). We have replaced this by the constant factor of 1/2. This is not 
quite proper since the period of the oscillation is a function of uj (and not just a constant). 
Near the cut-off frequency this error becomes the most apparent because the rest of the 
function in (ffOD is also a rapidly varying function (it is tending to a square root singularity). 
In the inset of Fig. we have plotted the ratio 
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FIG. 5. A plot of both linib^o 5'r(^/2.^) in ®, and 5^^(71/2,0;) in (H). Good agreement is 
obtained, but there is some discrepancy as shown in the inset, which shows the ratio of the two 
functions (0). This should be a constant (equal to one) if the two results are the same. The source 
of this discrepancy is discussed in the text. 

limfe_.ogf(7r/2,a7) ^ e-/3|^|/2j^(2A/3) 

5o^^(7r/2,a;) cosh (i/J^^p^ cot g) ' ^ ' 

as a function of uj for the same values of A, /3, and q used in previous plots. Near u; = the 
agreement is best (the ratio is near unity). But as oj — the agreement becomes progres- 
sively worse. Nevertheless, this rather crude treatment achieves reasonable agreement — only 
about 10% error at its worst. The agreement improves if we choose q 7^ 7r/2, but worsens 
as the product A/5 grows. (But A/? ^ 1 is the regime of interest.) 

This marks the end of the theoretical development. In the following section, we shall 
concentrate on various materials we believe are good candidates for observing the WZL. 
Specifically, we shall see that BO can exist in certain ferromagnetic Ising-like salts, with 
frequencies on the order of 150 GHz. 

IV. CANDIDATE MATERIALS 

We have identified four candidate materials for observing BO and the wzL in purely 
magnetic systems. The materials are all Ising-like FMs, and consist of chains of mag- 
netic ions, with effective spin-1/2, separated by spacer material. We focus mainly on 
C0CI2 -21120 (Ref. |12D, but give also a brief discussion on the potentially more promising, 
but less well characterized, CoCl2-2NC5H5 (Ref. 0), [(CH3)3NH]CoCl3-2H20 (Ref. §3), and 
[(CH3)3NH]FeCl3-2H20 (Ref. ||). 

A. CoCl2-2H20 

In C0CI2 -21120, the magnetic Co ions form chains along the c-axis. The coupling is 
ferromagnetic between ions in the same chain (we consider interchain exchange below). The 
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Parameter Value 
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9 


6.81 


5.49 


6.54 


7.49 




17.2 K 


3.17K 


4.14 K 


3.12K 


TABLE I. 


Relevant parameters 


for the candidate materials discussed in the text. 


A dashed 



entry ( — ) means that no value was given in the references. 

^CoCl2-2H20 (Ref. ||). 

bCoCl2-2NC5H5 (Ref. m). 

= [(CH3)3NH]CoCl3-2H20 (Ref. ||). 

'i[(CH3)3NH]FeCl3-2H20 (Ref. H). 



exchange anisotropy is such that the 6-axis is an easy-axis. The work of Ref. |T2| confirms 
unambiguously that the Ising-like spin-1/2 Hamiltonian of (§) describes this system very 
well. In Table |, we list the material parameters of this ferromagnetic salt. (We have taken 
the crystal 6-axis to coincide with the 2;-axis of the Cartesian coordinate frame.) We also 
list two antiferromagnetic interchain couplings. (We follow Ref. |T2| in neglecting the small 



non-Ising part of the interchain exchange.) We shall consider this interchain coupling in 
a mean field treatment by considering the total field at a given site to be the sum of the 
externally applied field and some internal field due to interchain exchange. 

Let us first neglect the interchain interaction and consider just a single chain in the 
presence of a static and homogeneous field along the z-axis. Then, the one-soliton approx- 
imation should be valid if the external field V = q^bH^^^ is less than the Ising exchange 
coupling . For the parameters given in Table |, this puts a restriction on the field strength 
of < 4T. For example, if we apply a 0.81 Tesla field — comfortably below this upper 
bound — then the Bloch amplitude and frequency are given from (^) as 

As = a^l^ ^ a, ^ = ^J^^ ^ 154 GHz. (75) 

This amplitude is small enough to impede the destructive influence of any scattering events, 
and the frequency falls within the capabilities of neutron scattering. This is therefore an 
encouraging result. 

Due to the antiferromagnetic interchain couplings Jf and J|, the material undergoes a 
three-dimensional ordering transition at about 17 K. Below this temperature, and in zero 
external field, there still exists ferromagnetic order within each chain, but the chains are 
ordered antiferromagnetically with respect to each other. As the field is turned on there are 
successive transitions from antiferromagnetic, to ferrimagnetic, and finally to ferromagnetic 
order at fields Hd and Hc2 respectively. At all times, the intrachain order is ferromagnetic. 
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FIG. 6. The three phases of C0CI2 -21120 for T < 17.2 K. The spin chains run perpendicular 
to the page and Jf and J| are the interchain antiferromagnetic couplings. These couplings can be 
considered internal fields, and so they too affect the Bloch oscillations. In fact, Bloch oscillations 



can potentially exist in all three phases. This figure has been adapted from Ref. 12 



This is depicted in Fig. ^, where we also list the interchain coupling values^ and the critical 
fields.EjtJ When determining the Bloch amplitude and frequency, one should also include 
these internal fields. For example, a Bloch amplitude of one lattice constant, which results 
from a total field of Htot = 0.81 T, can be realized in all three phases when these internal 
fields are taken into account. Explicitly, we have (for i^tot = 0.81 T) 

AFM : Hf^ = -i/""* - Ul/g^iB + 2J^/gi2B =^ i^'"* = 2.8 T, (76a) 



FIM : Hi"' = - U^/gfiB - 2J^/gfiB =^ H^""' = 3.6 T, (76b) 

FM : H\°' = i/""* + Ul/g^iB + 2J^/giJB =^ H^"' = 5.2 T. (76c) 

The notation -f^|°*, for example, denotes the total field at a chain with spin down, where 
"down" is defined as being opposite to the external field. Thus, in the ferromagnetic phase, 
all chains are spin up. In Fig. |^ we plot the resulting predictions for the Bloch frequency 
and (inverse) Bloch amplitude as a function of external field in all three phases. The antifer- 
romagnetic and ferrimagnetic phases show two curves because these phases have both spin 
up and spin down chains, and these chains each feel a different field. The dotted horizontal 
lines are bounds, outside of which the results become equivocal; near the upper bound, the 
total field becomes comparable to the Ising exchange constant J^; near the lower bound, 
the amplitude becomes too large, so that scattering effects should probably be taken into 
account. There is however, a fairly large intermediate range of over 400 GHz where the effect 
should be noticeable. We hope the present work can motivate some experimental investiga- 
tion into this material which seems to have disappeared from the attention of contemporary 
research. 
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FIG. 7. A plot of Bloch frequency and inverse Bloch amplitude in CoCl2-2H20 as a function of 
external field below the three dimensional ordering temperature of about 17 K. The discontinuous 
jumps in the curves are a result of transitions from anti-, to ferri-, and finally to ferromagnetic order 
of the chains relative to each other. Ferromagnetic order is always maintained within each chain. 
The horizontal dotted lines denote upper and lower bounds, beyond which the present analysis 
should not be expected to hold. 

B. Other Materials 

Another material we wish to mention here is CoCl2-2NC5H5 (Ref. ^Tj). This material 
differs from the one above only by the spacer material — pyridine molecules rather than 
water. Pyridine is a larger molecule than water and so the magnetic chains are further 
apart by about a factor of 1.7 (9.4A versus 5.5A for the water spacer). This material is 
thus a better one- dimensional material than CoCl2-2H20. The three dimensional ordering 
temperature is about 5.4 times smaller than it is in CoCl2-2H20 (3.17K rather than 17 K). 
The material parameters are summarized in Table |. Because the experimental work on this 
material seems less extensive than that on CoCl2-2H20, we have not made any predictions 
for BO in this material. But due to the reduced three-dimensional ordering temperature, 
this material should be a better candidate for observing BO and the WZL. Even if there 
really is no transverse anisotropy {Jy = Jx), BO could still be induced by simply applying a 
transverse field in addition to the one along the Ising-axis (see Sec. [II B 1| ). 

Most of the preceding paragraph applies even more emphatically to the final two materials 
listed in Table |, particularly to FeTAcH The small magnitude of the interchain couplings 
and the lower three-dimensional ordering temperature indicate that these materials may be 
quite suitable for BO. We have again chosen not to provide predictions for this compound 
since we believe further material characterization is necessary. 

In summary, we have shown in this section that there are a number of materials which 
may exhibit a dispersive soliton mode as well as BO. We have not discussed any AFM 
chains because we have been unable to identify any with the appropriate material param- 
eters such that the Bloch frequency and amplitude fall within experimentally accessible 
regimes. Should any such chains exist. Sec. |11 B 2| shows that BO may exist under an applied 
inhomogeneous magnetic field. 
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V. SUMMARY AND OUTLOOK 



In this work, we have shown that BO of magnetic sohtons occur in anisotropic spin- 
1/2 chains. Although we have focused primarily on biaxial Ising-like FMs, we have shown 
that BO can also occur in uniaxial FMs by applying a transverse field in addition to the 
longitudinal field. We have also shown that BO may exist in Ising-like afms by applying an 
inhomogeneous field. 

We have been mainly envisioning a neutron scattering experiment in this work because 
the dynamical structure factor shows clear evidence of the WZL; it contains sharp peaks at 
integer multiples of the Bloch frequency. At zero wave vector, all the peaks vanish except for 
the one at zero frequency and those at the Bloch frequency ±003- Thus, a measurement of the 
magnetization autocorrelation function (magnetic susceptibility) about the Bloch frequency 
should also detect the wzL. 

Several materials are promising candidates for observing BO and the wzL. We have cho- 
sen to focus our estimates on the one-dimensional salt C0CI2 -21120. Although this material 
is not an ideal one-dimensional substance (better ones have been identified above), BO of 
amplitude one lattice constant (about 3.6 A) and frequency of about 154 GHz are possible 
with applied fields of a few Tesla. The other materials we have mentioned are less well 
characterized than the one just described. However, the much smaller interchain coupling 
indicates that they are better one- dimensional samples — a statement further supported by 
their three-dimensional ordering temperature, which is much lower than CoCl2-2H20. It is 
possible that these materials contain only uniaxial anisotropy. If so, BO can be achieved by 
tilting the external field away from the Ising axis. 

A question we are currently investigating is what effect soliton-soliton interactions have 
on the WZL, as well as the related question on the influence of the higher-soliton states. The 
WZL and BO should survive if the number of solitons is conserved. In Ref. |T9|, a study is 
presented of soliton dynamics in the two-soliton sector (but in zero applied field). If one 
enforces periodic boundary conditions, then only even numbers of solitons can exist. But 
quantities calculated in the thermodynamic limit should be independent of the boundary 
conditions employed. Therefore, as expected, this work found practically the same result for 
the dynamical structure factor as Villain did working in the one-soliton sector. However, the 
two-soliton sector brings with it an opportunity to directly detect the coherent oscillation 
of solitons. For example, two solitons can form a bound state which should be identified as 
a magnon in spin-1/2 chains. Multiple magnon bound states can then be formed in which 
a cluster of adjacent spins are all flipped relative to the majority of the ferro magnetically 
aligned spins in the chain. Indeed, these are precisely the excitations measured in the 
work of Ref. O which concerns the optical excitation of multiple-magnon bound states. An 
enticing scenario exists if the ends of these clusters also undergo Bloch oscillation. Rather 
than having these excitations thermally created, as we have been assuming above, one can 
then optically create these excitations coherently by infrared radiation. The resulting Bloch 
oscillations will then also be coherent, and this may be detected, for example, by looking for 
coherent emission of magnetic dipole radiation in the microwave regime. In this scenario, the 
magnetic Bloch oscillator is an emitter of coherent microwave radiation. This is essentially 
the analog of the electronic BO experiments, were the charge carriers are optically excited 
with visible light, and the electron dipole oscillations radiate in the submillimeter regime. 
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There are no exciton effects in our spin cliains and so tlie detection of the magnetic radiation 
would be a clear signal of magnetic Bloch oscillations. This intriguing problem will be the 
subject of a future publication. 
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